ASYMPTOTICS OF MATRIX INTEGRALS AND TENSOR INVARIANTS 

OF COMPACT LIE GROUPS 



MICHAEL STOLZ AND TATSUYA TATE 

Abstract. In this paper we give an asymptotic formula for a matrix integral which plays 
a crucial role in the approach of Diaconis et al. to random matrix eigenvalues. The choice 
of parameter for the asymptotic analysis is motivated by an invariant theoretic interpreta- 
tion of this type of integral. For arbitrary regular irreducible representations of arbitrary 
connected semisimple compact Lie groups, we obtain an asymptotic formula for the trace 
of permutation operators on the space of tensor invariants, thus extending a result of Biane 
on the dimension of these spaces. 



1. Introduction 

Let G be a compact connected Lie group, and let (V\, p\) be an irreducible representation 
of G with highest weight A. Consider the matrix integral 



(TT Px (g)r(Trp x (g')r ■ ■ ■ (Trp A (/))^(Trp A (^)) 61 . . . (Trp\(g r )) dco G (g), (V 



G 



where {aj} r j =1 , {bj}j =1 are fixed sequences of non- negative integers and u G denotes nor- 
malized Haar measure on G. In the case that (V\, p\) is the standard representation of the 
unitary group G = U n , the integral Q is nothing else than the (ai, . . . ,a r ,b%, . . . , 6 r )-moment 
of the random vector 



(Tr(^),...,Tr(/),Tr(^),...,Tr(^)), (2) 

where g is chosen from U n according to Haar measure. It has been proven by Diaconis 
and Shahshahani in |DS] (see also |DEj ) that for n large enough, ([1]) coincides with the 
(di, . . . , a r , b%, . . . , 6 r )-moment of a vector of independent complex Gaussian random vari- 
ables of suitable variances, and consequently, the vector (j2]) converges in distribution to this 
Gaussian vector as the matrix size n tends to infinity. It has been observed ([BR], [St] ) that 
the Diaconis-Shahshahani result is based on the fact that the integral can be expressed 
as 

/ ]jTi(p x (gypTT(p* x (gyp du G {g) = Tr [{a ka ® a kb ){s,t)\ [v m a ^ rH] a) . (3) 

Here we write k a = YJj=ij a 3 an d define kf, analogously. au a denotes the obvious represen- 
tation of the symmetric group & ka on V® ka , and cr^ is its analogue on (V£)® kb 1 where the 
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representation (V x *,p* x ) is contragredient to (V\,p\). s G &k a has cycle type (l ai . . .r ar ) and 
t G &k b has cycle type . . . r br ). 

[Vx ka ® (V x f kb } G = {Te V® ka <g> (VZ)® k » : (pf ka ® (p*)® fcb )(^)T = T for all ^ G 67} 



is the space of invariants of the G-action pf a <g> (p* x )® kb . Diaconis and Shahshahani study 
the integral ([1]) for a sequence (G n ) of classical groups of increasing rank, fixing the param- 
eters a = (oi, . . . , a r ) and b = (bx, . . . ,b r ). Since they consider standard representations of 
classical groups, substantial information about the right-hand side of ([3]) is available since 
the epoch-making work of Weyl ( |Wy| ) and can be used to evaluate the left-hand side. 

In the present paper, we depart from the Diaconis-Shahshahani framework in two ways. 
Firstly, we consider arbitrary regular irreducible representations of arbitrary compact semisim- 
ple Lie groups. In this generality, one has a priori only poor control of the right-hand side 
of ([3]). So we will directly attack the left-hand side of this equation and thus obtain some 
asymptotic information about the spaces of invariants on the right-hand side. Secondly, we 
fix a group G and a representation V\, and let the parameters k a , kb of the tensor powers (and 
hence the moment parameters a, b) tend to infinity. This may be thought of as a thermo- 
dynamic limit of a particle system, rather than a random matrix limit. For our asymptotic 
analysis we will use techniques which were developed by Biane ([B]), Klyachko and Kurtaran 
QKK]), and Tate and Zelditch ( |TZj ). Our results should be compared to those of Biane (}B]) 
and Kuperberg QKJ). An asymptotic result of a different kind for growing tensor powers of 
a fixed representation has recently been obtained by Collins and Sniady in \CS\ Thm. 17]. 

Here is the set-up for our main theorems: Assume that the compact connected Lie group G 
is semisimple, and that the highest weight A of the fixed irreducible representation (V\, p\) is 
regular, i.e., is in the interior of a Weyl chamber. Fix a maximal torus T in G. Denote by W 
the Weyl group. Write t for the Lie algebra of T, and t* for its dual space. I := ker exp C t 
is the integral lattice, its dual I* : = {(pet*: <p(I) C Z} is the weight lattice. Let A* C t* 
be the abelian group generated by the roots, i.e. the root lattice, and write A := (A*)* for 
its dual. It is well-known that A* C /*, hence I C A, and that the group 11(G) := A/I is a 
finite abelian group. It can be regarded as a subgroup of T = t/I. Write tt for the canonical 
projection of T onto T /11(G). For any pel* write m^(/x) for its multiplicity in V\. Then, 
the set of all weights of V\ is M x := {pel*: m(/i) ^ 0}. 

Fix a smooth function / > on G and sequences a = (a%, . . . , a r ) and (3 = (Pi, ... , (3 r ) 
of nonnegative integers. For a positive integer set dj := cij(N) := Natj, bj := Nf3j (j = 
1, . . . , r) and write a = (a\, . . . ,a r ),b = (b\, . . . , b r ). Furthermore, we set 



a\ ;= y]aj, k a :=^2jaj, l a := }J 2 aj (4) 
j=i j=i j=i 
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and define |/3|, kp, lp, k a etc. analogously. Throughout the paper we will assume that k a = kp. 
We will study the integrals 



'TV 



W,<*) := / l[(Tr(p x (gi))f a if(g)dLU G (g) and (5) 

K N = K N (f,a,P) := / n(Tr(p A (^)))^ (Tr(p A (^)))^ f(g)du G (g). (6) 

- /G 3=1 

Writing $ + for the set of positive roots of (G, T), we define a polynomial k on t by 

k(x) = JJ (a,x), (7) 

Finally, we define A A £ Hom c (t, t*) 2* t* ® t* by 

Ax = dA^V ^ m M^®H- (8) 

Since G is assumed to be semisimple, and A to be regular, A\ is known to be positive definite 
(see [TZj ). 

Now we are in a position to state our main results, to be proven in Sections [3] - [5] below. 
Theorem 1. Assume that gcd{j : aj ^ 0} = 1. Then 

/Af(/ ' a) ~ ( 27 r/ Q iV)(^^/VdetA A 1 ^ { ) ' 

where d is the number of positive roots, p is half the sum of the positive roots, f G is the class 
function f G (g) = J f(x~ 1 gx)duG(x), and i>Nk a \ is the character on 11(G) determined by the 
weight Nk a \: 

v Nka x{h) = e 2 ^ lNk ^\ (10) 
•0^ £ A being a coset representative of h £ A/I = U(G). 

Theorem 2. Suppose that gcd{j : a 3 • ^ or (3j ^ 0} = 1. Then 

Remark 1.1. Note that in Theorem [T] the leading term vanishes if Nk a X = k a X is not 
contained in the lattice A* and if / = 1. In fact, in this case one has [V^® fca ] G = {o}, hence 
In(1,Q[) = by This is because the existence of a nontrivial invariant implies that 
there is a sequence pi, . . . , pt~ a of weights for the irreducible representation V\ such that 
Pi + ■ ■ ■ + Pk a = 0. Recall the well-known fact (see [TZ] ) that, if the highest weight A is 
regular, then the root lattice A* is spanned by the differences p — p! between two weights p, 
p! for V\. This implies that 

for any ip in the dual A of the root lattice A*, and hence k a X £ A*. 
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Remark 1.2. In view of ([3]), Theorems [I] and [2] give asymptotic formulae for the trace of 
permutations on the space of tensor invariants. Specializing to the identity permutation, one 
obtains the asymptotics of the dimension of these spaces. Specifically, taking in Theorem [1] 
atj = for j > 2 and ai = 1, / = 1, and assuming that A is in the root lattice A*, we obtain 
an asymptotic formula for the dimension of the space of tensor invariants, namely 

dimlVrf = |n(G)l(dimrA)jV/t(A " lp) (1 + 0(N'^)) 

L A J (27T)(^)/2iV(dim G )/2 v ^tA^ + ty ^ V >>■ 

This formula has been obtained by Biane in [B] (see also [TZJ). 
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2. Matrix integrals and tensor invariants 

Before turning to the proof of Theorems [1] and [2], we provide a short and self-contained 
proof of equation Q, on which the invariant theoretic interpretation of the integrals In and 
K N is based. Write H := [Vf ka <g> (V^) 9kb ] G . Then, the orthogonal projection onto H is 
given by 

tt h (T) = [ (pf ka ® (pir kb )(9)T du G (g), T G V® k ° ® (12) 
Jg 

For A G End G (V® ka ® (V^) m "), then, 

Tr(A\ H ) = Tt(Atth) = Tt(tt h A). 

We thus obtain 

Lemma 2.1. For any A G End G (Vf fca ® (V£)® kb ), we have 

Tr(A\ H )= [ Tr ((pf« ® (p* x r k *)(g) A) du G (g) = [ Tr (A (pf«®(p* x r k »)(g))dLU G (g). 
Jg Jg 

In view of this, ([3]) is implied by the following well-known lemma (see |DSj . [Raj): 

Lemma 2.2. Let V be a d- dimensional complex vector space, B G Endc(l / ), and s G &k 
a permutation of type (l ai 2 a2 . . . r ar ), hence k = k a = Y^j=\3 a i- Let &k act on V® k via 
Vk{s){® k =1 Vi) := ® k i=1 v ls -i. Then the trace of B® k a k (s) G End c (V r0fc ) is 

r 

Tr (B® k a k (s)) = IjTr(B J >. 

3=1 

Proof. Fix an inner product ( , ) on V . This induces an inner product on V® k , which is 
also denoted by ( , ). Let ei, . . . , be an orthonormal basis for V. Write T for the set 
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of maps from {1, . . . , k} to {1, . . . , d}, JF S for the restrictions to a subset S of {1, . . . , k}, 



ip ■ •<y v 



^e^f). Then {e v : ip G J 7 } is an orthonormal basis of V® k . Write 



r a j 



mi 

j=i i=i 



where {(f : z = 1, . . . , a?} are the cycles of length j in s. Furthermore, for any t G write 
[t] := {u — 1, . . . , k : z^£ ^ Then 

Tr(B® fc o- fc (s)) = {B® k ° k {s)e^ e,) = J2U (Be jt -i v , e jip ) 
tp<^T ipe^F j=i 

r a j r <%' 

= EnnnK)-v^)=nnE n^^-v^ 

j=l i=l j e [^j j=l i=l <fi£F [c j ] le[Cl] 

r a j r 

= nn Ti / > ] n Tii / > r 

3=1 i=l j=l 

□ 

3. Phase function for the matrix integral 

We start by rewriting the integral (jSJ) using Weyl's integration formula, assuming for 
simplicity that / is a class function: 



W,a) = ^J T UM P x^)) Naj f(t)\A(t)\ 2 dt, 



(13) 



where T C G is a maximal torus, di is Haar measure on T, normalized as a probability mea- 
sure and A(t) is the Weyl denominator. We define the following function on the complexified 
Lie algebra t : 

k(w) = m x (fj,)e 27r ^' w) , (14) 

where w = r + y/—lip G t© \/— It = t c , and linear forms in t* are extended complex linearly 
to t c . Note that the restriction of k to \/^Tt is essentially the character of V\. 

Let dip denote Lebesgue measure on t, normalized so that the fundamental domain T of 
the integral lattice I has volume 1. Then the integral Jjv can be written in the form: 



I N = -±- [ F(^lv) N f(v)\Mv)\ 2 d V , F(V=l<p) = l[k(V=ij<p) a > t (15) 
I JT j=1 



where the class function / is, through the exponential map, regarded as a function on t. 
Lemma 3.1. We have the inequality 

\F(V=l<p)\ < M0) H = (diml^) 1 " 1 . (16) 
Equality holds if and only if tp is in the dual lattice A of the root lattice A* . 
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Proof. Similar to the proof of Lemma 1.4 in [TZ]. Note that the assumption about the gcd 
is used here. □ 

By Lemma [3.11 the integral ffl~3l or ffl5l) is localized on ker7r = U(G) = A/ 1. Let g be a 
smooth cut-off function on T around the unit such that the support of g does not contain 
any element in ker7r other than the unit. Then, the translate Qh(t) = g(h~H) is a cut-off 
function around h G ker7r. Let (p G A. Then it is easy to see that 

k(V^ljip) = e 2 ^ j{x '^k(0), (17) 

which is not zero. Thus, around each (p G A, we can take a branch of the logarithm to define 
the following function $: 

r 

$(w) := y~] atj log k(jw), w = r + \^ltp G t c , (18) 

3=1 



where tp varies in a neighborhood of a point in A. Then, by Lemma 13.11 we can write the 
integral In as follows: 

I N = ±.e»*® J eN{ * {V=Iip) ~ m) 9h((p)f(<p)\A(cp) | 2 dtp (19) 

/iSker 7r 

plus a term of order 0(e _ciV ) for some c > 0. To compute each of the integrals in the sum 
in (fl9l) . we note that 

Afo + fa) = e 2n ^ T{p ^ h) A(ip) 
for each h G ker7r = A/I, where ^ 6 A satisfies exp(?/>/i) = h and p is half the sum of the 
positive roots. Furthermore, by ( fTTj) . we have 

HV^T& +ip h )) = HV^Tcp) + 2nV-l k a (\,?p h ). 
Therefore, we obtain: 
Lemma 3.2. We have 

I " = W\H v ^ h ) J e N *(^g(tp)f h (tp)\A(tp)\ 2 dtp (20) 

plus a term of order 0(e~ cN ), where g(ip) is a cut-off function around tp = 0, and fh{^p) = 
f(ip + iph) with a representative iph G A for h G ker n. 

Next, we compute the first and second derivatives of the phase function $ at points in 
ker 7r. 

Lemma 3.3. Any <p G A is a critical point of $. Furthermore, the negative of the Hessian 
o/$ ; H((p) : t — ► t* ; is given by H(<p) = (27r) 2 l a A\ which is independent of <p> G A and is 
positive definite, where l a and Ax are defined in (jl]) and (jHJ), respectively. 

Proof. That each ip G A is a critical point of $ is proven by the fact that 

rnx(jj)fji = 

because the left-hand side of the above is a jy-invariant vector and G is assumed to be 
semi-simple. The fact that the linear map A\ : t — * t* is positive definite is proven in 
[TZ] . □ 
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Therefore, what we need to do is to find asymptotics of each of the integrals in ( 120|) . each of 
which is an integral of functions supported around the origin. However, there is a difficulty: 
for each integral, the origin is the unique critical point for the complex phase function $, 
but the Weyl denominator A vanishes at the origin. In the next section, we will use ideas of 
Biane ([B]) and Klyachko-Kurtaran ( [KK ] ) to circumvent this problem. 



4. The method of Biane and Klyachko-Kurtaran 

In this section we study the following integral, which is a slight generalization of the 
integral on the right-hand side of (|20|) : 

J N := [ e N ^g(t)\A(t)\ 2 dt, (21) 



JT 

where dt is normalized Haar measure on the maximal torus T and A denotes the Weyl 
denominator, g is a compactly supported smooth function which does not vanish in the 
unit element of T. We fix a jy-invariant inner product on t such that the volume of the 
parallelotope determined by an orthonormal basis equals 1. We assume that the complex 
valued smooth phase function $, which we view as a function on the complexified Lie algebra 
t c , satisfies the following conditions: 

(i) Re($(v / ^Tv 5 )) < $(0), with equality if and only if (p = 0; 

(ii) the origin is a critical point of $; 

(iii) H = — <9 2 $(0) : t —>■ t*, the negative of the Hessian of $ at the origin, is positive 
definite; 

(iv) in the Taylor expansion 

$( V ^1 V9 ) _ $( ) = -(H<p,<p)/2-V=i &(<p) + RM (22) 

up to fourth order, 6 is a real valued homogenous polynomial of degree 3 

(v) the linear map H commutes with the action of the Weyl group W. 

Note that the phase functions of the integrals 1^ and K N (see Sections [3] and E]) satisfy these 
conditions. We are now in a position to state the main result of this section: 

Theorem 4.1. Under the above conditions on $ and g one has 

, 2 x(dimG)/2 (2 n )d a (n) e NZ(0)\ W \ 

Jn = g j (2.) fl (0)e_ \W\ <H _ lp) (1 + 0{N _ m)l p3) 

where \ W\ is the order of the Weyl group W, d is the number of positive roots, p is half the 
sum of the positive roots, and the polynomial k on t is defined in (|7j) . 

Proof. First of all, we normalize Lebesgue measure on t so that the volume of the fundamental 
domain T of the lattice I equals 1, and write 

J N = e N ^ [ e N W^-*Wg&)\A&)\ 2 dp. (24) 

We may regard g as a function on t with arbitrarily small compact support around the 
origin, since by (i) the integrand in (1241 is bounded by e~ cN (with a constant c > 0) outside 
a compact neighborhood of the origin, (i) and (iv) imply that for fixed a > one can choose 
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b > such that Re$( v /T T^) - $(0) < -b(H(p, <p) for \<p\ < a. Substituting ([22]) into J N 
and changing the variable ip to N~ l l 2 ip, we have 

J N =A MrkG)/2 e ^(0) 



As in [B], |TZ] . using the identity A = riae<i> + (e 7r ^^- a — e 7r, ^i Q: ) ) it is easy to see that 
A(N~ 1/2 (p) = (2irV^T) d N- d/2 K(<p)(l + 0(N~ l \ip\ 2 )). 

We note that g(N~" l l 2 <p) = g(0)(l + 0(N-^ 2 \ip\)). Substituting these formulas into ([25]) and 
introducing a smooth cut-off function \ sucn that X — 1 around the support of g, we obtain 

_ (2n) 2d g(0)e N ^ 

f +(rkG)/2 _ (26) 

x / ^(^.i^^-v^M+^r./v)^^-!^^^^^ (l + otiv- 1 / 2 )). 



Here, we note that g and its derivatives are bounded on t, and that the exponential in the 
integrand is bounded by e~ b ^ H<p,<p ^ if \tp\/N l l 2 < a. As indicated above, we may assume that 
g(N~ 1 ^ 2 (f) = for \ip\/N 1 ' 2 > a. 

Next, as in |TZ], we divide the integral in (1261) into several parts as follows. We set 

E N (ip) ■= e -(H^)/2-^iN-^e( v ) and write 

where the integrals Ij(N), j = 1, 2, 3, 4 are given by 



h(N) = J E n ( V )\k( V )\ 2 dip, 
I 2 (N) = j ^M(e^- 1/ ^-l)x(iV- 1/ V)l^)| 2 ^, 
h{N)= j S iV ( V p)e iV ^^ 1/2 ^(l-x(iV- 1/ V))x(iV- 1/ V)l«(^| 2 ^, 
h{N) = [ E N {<p)(x{N-^ip) - lMip)\ 2 dip, 



(28) 



where we used the relation x(N 1 ^ 4 ip)x(N 1 ^ 2 ip) = x{N 1//4 y?) f° r sufficiently large N. 

Now, the integrand in I2(N) vanishes for \cp\ > cN 1 ^, and thus we have e NR ^ N 1/2</3 - ) = 
0(1) and NR^N- 1 / 2 ^) = |^| 4 0(l/iV). Hence, ( e NR ^ N ~ 1/2 ^ - 1) is bounded by c|y?| 4 0(^). 
But we have \En(^>) \ = e - '^'^" 2 , and hence \ip\ 4 EN(ip) is integrable uniformly in N. Thus 
we have I 2 (N) = 0(1/N). As to I 3 (N), the function E N (ip)e NR ^ N ~ 1/2 ^ is dominated by 
e -b{Htp,<p) wnereV er the integrand does not vanish. Since x{N~ l l A ip) = 1 for \cp\ < cN 1 ^ for 
some c > 0, we easily have I%{N) = 0(N Cl e~ C2Nl 2 ) for some ci,c 2 > 0. Similarly, we have 
I 4 (N) = 0(N c ^e- C2Nl/2 ). 



< Hx ^/ 2 \K( X )\ 2 dx 



ASYMPTOTICS OF MATRIX INTEGRALS 9 

Finally, we consider the integral h(N). Note that e - ^ 1 @ ^) N ' 1 ' 2 = 1 + O^ipf- /N 1 ' 2 ) 
since Q(<p) is real. Thus, invoking the identity 

VdetiJ 

known as Mehta's conjecture and proven in |Mcj and |Up] , we have 

h (N) = ! e-^^\K{ V )\ 2 d^l + 0{N-" 2 )) = ( 2 K) irkG)/2 m< H^ti {1 +0(iNr i/2 ))> 
J V det 

(29) 

which completes the proof. □ 

5. Proof of the main results 

Theorem [1] is a direct application of Theorem 14.11 with Lemma 13.31 to the integral on the 
right-hand side of ( !20l) . To prove Theorem [2} we proceed as in Section [3] and use the Weyl 
integration formula to obtain 

K N = ^J T J(tff(t)\A(t)\ 2 dt, J(t) = nTr(pA(^)rTr(p A (^)f . (30) 

We also apply Theorem 14. II to find asymptotics of the integral K N as follows. As in Lemma 
13.11 we have 

\J(V-l<p)\ < M0) H+I/3 ' = (dimV A ) H+l/31 , (31) 
with equality if and only if ip G A. We define 

s 

*W = Y)- a J log k ^ w ) + Pi log k ^ W ^ w = r + y/=T<p G t c , (32) 
i=i 

around each ip G A. Next, we need to compute the Hessian of \& at <£> G A. We note that, 
for ip G A, we have 

fc(V=I?» = e 2 ^< A ^>A;(0), (dk)(V=lj<p) = 0, 
9 2 (logA;(v / ^I^)) = -47r 2 j 2 A A . 
Therefore, we have the following lemma: 

Lemma 5.1. Each (p G A is a critical point of^>. The negative of the Hessian of^, denoted 
D(<p), is given by D(<p) = (27r) 2 (l a + l(3)A\, and hence D = D(<p) does not depend on ip G A, 
and is positive definite. 

This time, the term involving 27i\/^lj{\,iph) disappears because of the assumption 
jotj = Y^JPj- We thus have 

Kn = W\H Je N ^^g(<p)f h (<p)\A(<p)\ 2 d<p, 

where fh is defined in Lemma 13.21 The assumption of Theorem 14.11 is satisfied by the phase 
function and applying it, we obtain Theorem [2j 
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